Energy momentum tensor, stability, and the D-term of Q-balls 
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We study the energy-momentum tensor of stable, meta-stable and unstable Q-balls in scalar field 
theories with U(l) symmetry. We calculate properties such as charge, mass, mean square radii 
and the constant d\ ("D-term") as functions of the phase space angular velocity u. We discuss the 
limits when uj approaches the boundaries of the region in which solutions exist, and derive analytical 
results for the quantities in these limits. The central result of this work is the rigorous proof that di 
is strictly negative for all finite energy solutions in the Q-ball system. We also show that for Q-balls 
stability is a sufficient, but not necessary, condition for di to be negative. 
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I. INTRODUCTION 

The matrix elements of energy momentum tensor T t 
(EMT) m contain basic information, such as mass 
and spin of the particle Q, and the constant d\. This 
constant denotes the value of the corresponding form fac- 
tor at zero-momentum transfer, and is the only experi- 
mentally unknown "charge" associated with the EMT. 
Its physical meaning was uncovered in Q where it was 
shown that d\ is related to the spatial distribution of in- 
ternal forces. This is analog to the interpretation of the 
electric form factor as the Fourier transform of the elec- 
tric charge distribution [BJ (and subject to the same type 
of limitations 0,0] due to relativistic corrections). 

The EMT form factors can in principle be studied 
through generalized parton distribution functions 0, Q 
accessible in hard exclusive reactions such as deeply vir- 
tual Compton scattering [§-[l2|]. In this work we shall 
loosely refer to d\ as the D-term, although they coincide 
strictly speaking only for asymptotically large renormal- 
ization scales [13[ . Theoretical studies of EMT form fac- 
tors were presented in chiral perturbation theory, lattice 
QCD, effective chiral field theories or models |13l - [23| . 

It is a striking observation that in all theoretical stud- 
ies d\ was found negative — for pions, nucleons, nuclei. 
Results from chiral soliton models [mjl9| suggest a con- 
nection of the sign of d\ and stability of the particle. The 
naturally emerging questions arc: could the negative sign 
of d\ be a model-independent feature, a theorem? And, 
what is the precise relation of d\ and stability? 

This work is devoted to the study of the EMT of Q- 
balls [13, HH with the aim to shed further light on the 
sign of di and its connection to stability. Q-balls are non- 
topological solitons in theories with global Abelian [25[ or 
non-Abelian |26j symmetries, and have been discussed in 
a variety of approaches with a wide range of applications 
in particle physics, cosmology, and astrophysics (27l - l45j . 

The Q-ball equations of motions admit stable, meta- 
stable and unstable solutions. This makes them an ideal 
ground for our purposes to explore the connection of the 
stability of an object and its D-term. In this work we will 



be interested in the ground state Q-ball solutions [25| . 
The equations of motions admit also radial excitations 
[37| which will be subject to a separate work [46] . 

The stability of Q-ball systems was studied in many 
works [37l-l4ll]. But this is, to the best of our knowledge, 
the first time this issue is addressed from the point of 
view of the EMT, and the D-term. In particular, we 
present the first rigorous proof in a dynamical system 
that d\ must be negative, and clarify the connection of 
d\ and stability. The outline of this work is as follows. 

In Sec. [nl after a brief review of Q-balls, wc derive the 
expressions for the energy density Xbo(r), pressure and 
shear forces, p(r) and s(r), related to the stress tensor 
Tik(r), and prove analytically that exact solutions of the 
Q-ball equations of motions satisfy the stability condition 
which is a consequence of the conservation of the EMT. 
We show that the stability condition is equivalent to the 
virial theorem, as it is the case also in soliton models of 
the nucleon 

In Sec. IIIII we study Q-ball properties such as charge, 
mass, mean square radii, and the D-term in a chosen 
potential as functions of the angular velocity uj in the 
U(l)-space in the region oj mnl < u < w max in which the 
equations of motion admit finite energy solutions. An 
interesting observation is that among the quantities we 
study d\ varies most strongly with to. 

In Sec. llVl andlVlwe then focus on the behavior of the 
Q-ball properties as uj approaches the boundaries of the 
region in which finite energy solutions exist, i.e. the limits 
uj — > oj m - nl and u> — > w max respectively. We derive in both 
cases analytical results which describe the behavior the 
different quantities in these limits and which are fully 
supported by the numerical results. 

In Sec. I VII we formulate two independent proofs that 
d\ of Q-balls is strictly negative, and clarify the relation 
to stability. We show that stability is a sufficient but not 
necessary condition for d\ to be negative. 

The Sec. I VIII contains a summary of our findings and 
the conclusions. Some technical details and supplemen- 
tary results are discussed in the Appendix. 
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II. Q-BALLS 



A. Stability criteria 



We consider the relativistic field theory [25[ of a com- 
plex scalar field defined by the Lagrangian 



C 



1 



(1) 



The potential V > is defined in terms of the positive 
constants A, B, C (with AAC > B 2 to guarantee V > 
for <f> == 0) as follows 

V = A ($*$)- J B($*$) 2 + C($*$) 3 . (2) 

The Lagrangian is invariant under global U(l) symmetry 
transformations $ — > <I>e IQ , $* — > <I>*e~ IQ; with a real. 
This system admits non-topological soliton solutions (2f| 
which in the soliton rest frame are given by 



$(i, x) — exp(iujt) 4>{r) , r = \x \ 



(3) 



The Euler-Lagrange equations of the theory in (TTJ) imply 
for the radial field 4>(r) the following differential equation 
(here and in the following primes denote differentiation 
with respect to the argument) 

// (r) + % , (r)+wV-V"(^)=O, (4) 
r 

which is subject to the boundary conditions 

0(0) = 0o = const, 0'(O) = 0, 

0(r) -> for r — y oo . (5) 

The Noether theorem applied to the global U(l) sym- 
metry implies the conserved charge 



Q 



d 3 x p c h(r) , Pch(r) = u) (t>(r) z 



(6) 



The sign of w determines the sign of the charge Q. In the 
following wc assume u> > without loss of generality The 
presence of a continuous global symmetry is essential for 
the existence of the soliton. More precisely, finite energy 
solutions exist for u in the range 



2 2 2 

J min < w < w max ' 



(7) 



with 



2V{<j>) 



= 24 1- 



B 2 \ 
4AC ) 



> 



< ax = V\4>) 



= 2A 



(8) 



i=0 



From (jHEJ) we obtain for 0(r) the following small- and 
largc-r behavior (the dots indicate subleading terms) 

4>{r) = O + (v'(M - u 2 <f>o) + ■ ■ ■ small r, (9) 



cxp 



+ . . . large r. (10) 



The constants 4>o and Coo are known, of course, only after 
solving the boundary value problem ((HIS]). 



Solutions for u satisfying the existence condition ([7]) 
can be classified as (a) stable, (b) meta-stable, and (c) 
unstable Q-balls, see e.g. [4l| for an overview. 

(a) If M denotes the mass of the soliton, and m the 
mass of the field $, which is m = s/2A = w max , then the 
absolute stability condition can be expressed as [29[ 

M < mQ , m = w max . (11) 

(b) Meta-stable solutions do not satisfy pT|) but are 
stable with respect to small fluctuations, and satisfy a 
weaker "classical stability condition" [24|, H^] which can 
be formulated in the equivalent ways 



M 
dco V Q 



dQ d 2 M 



(12) 



i.e. a critical ui c (extreme charge Q c ) exists at which the 
slope (curvature) of the quantities in (|12[) changes. 

(c) Solutions satisfying neither the stronger condition 
(jTTJ) nor the weaker condition p^|) are unstable. 



B. The EMT of Q-balls 

For the theory defined by the Lagrangian ([TJ the 
canonical energy momentum tensor 

T ^ = W^) d ^ + W^) d ^- 9 ^ c (13) 

is symmetric and static. The energy density, which de- 
fines the mass M = J d 3 xToo, is given by 

Too(r) = \ u 2 4>{r) 2 + \ 0'(r) 2 + V(4>) . (14) 

The Tot components vanish, i.e. the Q-ball has spin zero. 
(Of course, Q-ball solutions can be assigned a non-zero 
spin by means of appropriate projection techniques |47j.) 
Finally, the components describe the stress tensor 



Ti 



^-ls ij )s(r)+6 ij p(r) 



with the distribution of the shear forces, s(r), and pres- 
sure, p(r), given by 

s(r) = (f>'{r) 2 
1 

2" " v/ 6 
The dimcnsionless constant d\ is defined through the 
stress tensor Tij and can be expressed in terms of 
s(r) and p(r) (cf. Sec. IIIC|) as follows 



(16) 

= -0J 2 4>(r) 2 - l 4>'( r ) 2 -V^). (17) 



47T f°° 

di = - — M / dr r 4 s(r) , 
3 Jo 

/>oo 

= 57r M / dr r 4 p(r) . 
Jo 



(18) 
(19) 



The large-r asymptotics (|TU1) ensures that all integrals 
(which define M, Q, di, mean square radii, etc) exist. 
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C. Consequences from conservation of EMT 

For a static EMT d^T^ = is equivalent to V*7y = 
which, using the decomposition (fT5)). implies [i| 



(20) 



In order to prove that (|20|) holds for Q-balls, we insert 
s(r) and p(r) from Eqs. ([TBI [T7| into (f2T)|). which yields 



r 3 



(21) 



= 4>'{r) U"(r) + - 0'(r) + w 2 0(r) - V'{4>)j = 

due to the equations of motion in Eq. Q. 

From (|20p one can derive the equivalent representa- 
tions (fTSl IT91 for dn in terms of s(r) and p(r), and other 
general relations [Ta] . For instance, multiplying ([2"U|) by 
r and integrating (by parts) over r from zero to infinity 
yields the "stability condition" (also referred to as "von 
Laue-condition," see HH and references there in) 



dr r 2 p(r) = . 



(22) 



In order to prove (|2"2"j) for Q-balls we integrate by parts 
(primes denote derivatives with respect to the arguments, 
the finite upper integration limit R is for later purposes) 



(23) 



R 




R R 


J dr r 2 p{r) = 


y P{r) 


.7 










Next we notice that 



p'(r) 



--0"(r)+ W 2 0(r)-y'(0))0'(r) 

|0W'(r)--0'(r) 2 
3 r 



2 

3^3 



r^'(r)^ 



(24) 



where we used uj 2 <fi — V'(<fi) = —<f)"(r) — - </>'(r) in the 
first step which holds due to the equations of motion Q . 
Hence, using s(r) = 0'(r) 2 , we obtain 



/ 



dr r 2 p(r) 



t ( . . 2 . . 
y I P(H + 3 s(r) 



(25) 



The small- and large-r behavior of the solutions in © 
guarantees that the lower and (after taking R — > oo) 
upper integration limits in (|25[) vanish which proves (|22p . 

It is instructive to prove ([22]) independently as follows. 
Let 0(r) be a Q-ball solution with charge Q and mass 
M, which we rewrite by means of (0 ITU in terms of the 



"charge," "surface," and "potential energies" (where we 
leave the number of dimensions D — 3 general) 



1 1 

M = — + — E sm { + E pot 



(26) 



d D x(j)'(r) 2 
d D xc/)(r) 2 



E pot = d u xV{<t>) 



E ch = y- 



(27) 



We consider dilatational variations <fi(r) — > 4>(\r) of the 
solutions with a positive parameter A. Substituting in 
the integrals in ([27]) x — > Xx yields 

I _ , n 1 
2 



M(A) 



Esmf A 2 15 



-Epot A D 



(28) 



M(A) has a minimum at A = 1, because 4>(Xr) for A = 1 is 
a Q-ball solution which minimizes the energy functional. 
We obtain, using again ([5]) and the definitions in (1271) . 







1 dM(X) 
D 



OX 



1 



A=l 
2 



~\ Ech 



2- D 
2D 



E s 



surf 



E 



pot 



2D 



^{rf -V^) . (29) 



For D = 3 we identify the expression (|17p for p(r) in the 
curly brackets of (|2"9"|) which completes our alternative 
proof of (|2"2"|) . Eq. (|29p is known as virial theorem [30| . 
Notice that (|29p can be used to eliminate, for instance, 
the potential energy term from ([26|. leading to 



M 



1 

D 



surf 



(30) 



As a last application of (|2TJ|) wc integrate this equation 
over r from zero to infinity. This yields the relation [l|| 



p(0) = 2 / dr 



s(r) 



(31) 



which provides a helpful cross check for numerical cal- 
culations, and implies the following interesting relation: 
inserting in (|3"Tj) the expressions ((TBI [T71) yields 



i^ 2 -y(0 o ) = 2 



dr 



</>'(r) 2 



(32) 



o 



This is interesting, because the left-hand-side depends on 
</>o only while the right-hand-side is a functional of 4>'(r) 
where <fio drops out. Below in Sec. IIII Al we will discuss 
the physical interpretation of (|32|) . 



D. Relations among Q-ball properties 

Further interesting relations among different Q-ball 
properties follow from combining (|6l [T4J [TBJ [T7|) as 



7oo(r) + p(r) = u) p ch (r) + - s(r) 



(33) 
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Integrating (|5"3"|) over d 3 a; we recover (f3"0"l) for D = 3 (the 
derivations are equivalent, but (|30p elucidates the rela- 
tion of the factor ^ to the dimensionality of the space). 
Next, we define the "surface tension" 7 and the mean 
square radius (r 2 ) of the shear forces s(r) as follows 



III. GROUND STATE Q-BALLS 

In this Section we discuss the ground state properties 
of Q-balls in our potential (|2|) for different values of u>. 



00 J dr r 2 s(r) 

7 = J dr s(r) , (r 2 s ) = ^ 

J dr s(r) 



(34) 



= Jd 3 xs(r), Eq. (H7J), 
2 ^ 7 which is what one 



Thus, the surface energy E sur { 
can be written as E SUI f = Air (r 
expects for a spherical object with a well-defined surface 
and radius (r 2 ) 1 / 2 . In Sec. IIVI we will see that these 
notions make sense for Q-balls in a certain limit. 

Finally, we weight (|3"3"|) with r 2 , and integrate over d 3 :c. 
This allows to express d\ in terms of other properties as 



= ^[uQM(r 2 Q ) 



M 2 (r|; 



(35) 



with the mean square radii of energy and charge densities 
defined as 



(r 



Jd 3 x r 2 T 00 {r) 
Jd 3 xT 00 (r) 



Jd 3 xr 2 p ch (r) 
J d 3 xp ch {r) 



(36) 



E. Parameters and numerics 

In our numerical study we fix the parameters as 

A = 1.1, B = 2.0, C=1.0. (37) 

(for which in j37j radial Q-ball excitations were found; 
the latter originally motivated our study, but will be dis- 
cussed in a separate work (4(| ) . This yields the following 
range of allowed w-values 



0.2 < uj z < 2.2. 



(38) 



The parameter set (A, B, C) could be assigned physical 
units, say (GcV 2 , GeV°, GcV~ 2 ). Then u, M would be 
given in GeV, mean square radii in GeV -2 , etc. But for 
simplicity we will work with dimensionless quantities. 

The numerical method is as follows. For a given uj the 
differential equation (|4]) is solved with slightly shifted ini- 
tial conditions 4>{e) = <j) e and 4>'{e) = ^(V (4> e ) — Lo 2 4> e )e 
with numerical parameters e in the range 10~ 10 to 10~ 4 . 
We checked that the results do not depend on e. Finite 
energy solution are found using the shooting method by 
varying the initial value <p £ until 4>{r) — > at large r. 

The quality of the numerics is monitored by testing 
that (i) the differential equation ([20|) holds, (ii) the sta- 
bility condition (f2"2"j) is valid, (Hi) different expressions 
for di in (fT8l [T9l [35]) yield the same result, (iv) the same 
value for p(0) follows from (flTl l3Tj) . We find a relative 
numerical accuracy of O(10 ~ 6 ) or better. 



A. Effective potential U e ff and ef> 

Identifying r — > t and <p(r) — > x(t), the equation of 
motion Q can be read [25[ as the Newtonian equation 

(39) 



x(t) = F fric - VC/ eff (x) 

1 fllC 



1 



■ x(t) , U cS = -w 2 x 2 - V{x) . 

describing the motion of a particle of unit mass under 
the influence of the time- and velocity-dependent friction 
Ffric in the effective potential U e ff shown in Fig. [1] The 
initial and boundary values (|5|) mean that at t = the 
particle starts from the position xq with zero velocity, 
and comes to rest in the origin x = after infinite time. 
Thus x(t) > and the particle never stops at finite t. 
This implies decreasing monotony of the ground state 
fields, 4>{r) > and 4>'(r) < for < r < 00. 

The pressure (fT?) is given at r = by p(0) = C/ e ff(</ , o) 
and the condition ([7]) gua rantees the existence of a region 
of <f> with t/cff(</>) > [25[ . This proves that 



p(0) > 



(40) 



Now also the physical interpretation of (|3"2"T) is evident. 
The left hand side of (|32|) is the initial potential energy. 
The right hand side of (J22]) is the work W = J F flic dx 
the particle does to overcome the friction before coming 
to rest at x = with zero effective potential energy. 
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FIG. 1: The effective potentials C/ ff(c/>) = \^ 2 <t> 2 - V(<j)) as 
functions of <j> for selected values of u> 2 in the range ()38|) . The 
circles show the initial values c6o for each oj 2 , which lie on a 
curve starting and ending at the limiting values = 0.2 

and c^max = 2.2 (marked by open circles). The special values 



1.9 and lj 



2 

abs 



1.55 are discussed in text. 
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Solutions exist for all ui in the range uj min < uj < ui max . 
With our numerical method, we were able to find solu- 
tions in the subinterval 0.216 < ui 2 < 2.195. 

The effective potentials U e s(4>) = \uj 2 <j) 2 — V{<f>) are 
shown in Fig. Q] for selected uj including the limiting cases 
w m in and w max . On each of the E/ c ff(0)-curves in Fig.[T]we 
marked the initial conditions c/jq which solve the bound- 
ary value problem (jH [5]). The U e s(^>o) f° r different u 
lie on a curve which exhibits a global maximum close to 
w abs ~ 1-55, and changes the curvature around ui 2 ~ 1.9 
(these frequencies will be discussed in detail in Sec. lIII Cl) . 
This curve starts and ends at the limiting points 



lim 



00 =0. 



lim 




1 



(41) 



with the potential U e s{4>o) — > in both cases, see Add. [XI 
U c ff((f>o) as function of 4>q is n °t unique for 4>o >1. 



B. Solutions 4>(r) and densities 

In this section we describe the results for <j)(r) and the 
various densities. Some of our observations concerning 
the behavior of the densities in the limits uj — > w m in,max 
will be made more rigorous in Sees. ITVl andlVl 

The ground state solutions <f>(r), which are uniquely 
determined in terms of the initial values 4>q discussed in 
the previous section, are shown in Fig. [2^, as functions of 
r for selected values of u) in the range 0.216 < uj 2 < 2.195 
our numerics can handle. We have chosen a logarithmic 
r-scale to better show the features of all solutions in a 
single plot. On a logarithmic scale the <j)(r) are nearly 
constant for r < 0.1, and have the small- r behavior ©. 
Their large- r asymptotics agrees with (fTU| . 

With decreasing uj the solutions <p(r) remain nearly 
constant at their initial values in a region < r < Rq 
and form increasingly long plateaus from which they then 
drop down to their large-r asymptotics (|10[) over decreas- 
ingly narrow transition regions with thicknesses <C Rq- 
Here Rq can be understood as the "size" of the Q-ball, 
which will be defined below more accurately. In the limit 
u — > ujmin the field <p(r) — > <fio O(i?o — r ) where Rq — > oo 
and <f>Q -> 1, sec (|4"Tj) and App. This behavior ("thin- 
wall limit") can be strictly derived [25]. In the oppo- 
site limit, as uj increases, the solutions <f)(r) become more 
wide-spread and their magnitude decreases, see (|4i~j) and 

App.m 

Fig. [2b shows the charge densities p c h(r) = ui(j)(r) 2 
as functions of r. Also the charge densities exhibit for 
small uj extended plateaus in the region < r < Rq 
inside the Q-balls, and drop abruptly to zero outside. For 
uj — > w max the charge densities become more wide-spread 
and their magnitudes show an overall decrease. 

Fig. [2t shows the energy densities Too i r ) > which look 
qualitatively similar to charge densities for uj > 1. But 
for uj < 1 the energy densities start to develop a "bump" 
around Rq , and as uj approaches w m i n the bump becomes 
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FIG. 2: Field 4>{r), charge density p c h(r), energy density 
Too(r), pressure p(r), shear force distribution s(r) vs. r for 
Q-ball ground state solutions for selected values of uj. 
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a characteristic "spike." The reason for that is that for 
u) > 1 the Q-balls arc "diffuse" objects, while for uj < 
1 they start to develop a more and more well-defined 
"edge." In fact, asw-> w m in the notions of a "surface" 
and "surface tension" become better defined [25|. The 
characteristic bump/spike structure in Tqq(t) at r ~ i?o 
reflects the contribution of the surface energy. In the 
limit u> — > w ma we find Tbo(r) — > 0. In the limit w — >• 
w m in we have T o(V) — > const for < r < R with a 
surface energy contribution proportional to 5(r—Ro) with 
i?o — > oo. 

Fig. [2ji shows the pressures p(r) as functions of r. For 
all ground state solutions the pressures are positive "in- 
side" and "negative" outside in agreement with general 
expectations (|75|) . We are now in the position to provide 
an exact definition of the scale Rq as the point where 
the pressure vanishes, i.e. p(Ro) = with < Rq < oo. 
The stability condition is fulfilled because the integrals 
J Q R °drr 2 p(r) and j^dr r 2 p(r) have opposite signs and 
precisely cancel. Numerically the sum of these 2 con- 
tributions normalized with respect to the sum of their 
moduli is of 0(1O -6 ) or smaller. In the limit cj — > w max 
we find p(r) — > 0. For uj — > uj m i n we obtain p(r) — > const 
for < r < Rq with a surface energy contribution pro- 
portional to — 5(r — Rq) and a diverging Q-ball size Rq. 

Fig. [5^ shows the shear forces s(r) which are best suited 
to discuss the concepts of "diffuseness" or "edge" . From 
(HI [lOj) we see s(r) — > as r — > or r — > oo, and 
from Eq. (|I6|) we see it is an evidently positive quan- 
tity, i.e. s(r) must have a global maximum somewhere. 
To determine the position of this maximum consider 
s'(r) = 2<f>'(r)<f)"(r). Now due to the monotony property 
of <fi(r) discussed in Sec. IIII A[ we have 4>'{r) = only at 
r = and at infinity. I.e. the maximum of s(r) coincides 
with the point where 4>"(r) = 0. This change of curva- 
ture occurs in the vicinity of the "edge" of the Q-ball, 
and in the limit u> — > w m i n precisely at r = Rq where s(r) 
becomes proportional to a 5(r — Rq) with a coefficient 
related to the surface tension. As u> approaches w max the 
shear force distribution becomes wider and wider, which 
indicates a more and more diffuse "edge" of the Q-ball. 



C. u> c and oj a bs 

The frequencies uj c and w a bs were discussed in the se- 
quence of Eq. (fl"2")) , and mentioned in the context of 
Fig. UJ In this Section we discuss how they appear in 
the numerical results. At the frequency uj — lu c « 1.38 

1. (M/Q)(u>) has a global maximum, Fig. [3Ji, 

2. Q(lu) has a global minimum, Fig. [3J}, 

3. M{Q) has a branch point at Q c = Q(uj c ), Fig.[3J:. 

The frequency ui c and charge Q c define classical stabil- 
ity. For uj < uj c we have (M/Q)'(oj) > 0, and Q'(u) < 0. 
For Q > Q c we have M"(Q) < 0. These are equivalent 



M/Q (a) 




FIG. 3: (a) M/Q as function of u, with (M/Q)'(oj) > in 
the range ui < uj c ~ 1.38 (solid line), (b) Q as function of u), 
with Q 1 (uj) < for uj < uj c (solid line), (c) M vs. Q, with 
M"(Q) < in the branch denoted by the solid line. The solid 
(dotted) lines correspond to the region of classically stable 
(unstable) Q-balls, see Eq. (fT2j) . 
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FIG. 4: (M — mQ)/M as function of uj which exhibits a 
maximum at uj c . For uj < uu a bs we have M — mQ < and the 
Q-balls are absolutely stable. For uj a b B < uj < uj c the Q-balls 
are meta-stable, and for uj c < uj < c^max they are unstable. 
(M — mQ) approaches zero as uj — > oj m ax which is indicated 
by the vertical line. 



criteria for the stability of Q-balls against small fluctua- 
tions. In Fig. [3] the branches of classically stable Q-balls 
are shown as solid lines, while the unstable branches are 
depicted as dotted lines. 

Classical stability is a necessary but not sufficient con- 
dition for stability. For a Q-ball to be absolutely stable it 
is required M < mQ, sec Eq. ([TTj) . Fig.|4]shows M — mQ 
normalized with respect to M as function of uj. The 
quantity (M — mQ)/M is negative for uj < w a b s ~ 1.245 
(and has a global maximum at uj = uj c which follows from 
the fact that M(uj) and Q(w) have extrema there). Thus, 
for uj < w a bs the Q-balls are absolutely stable. In the re- 
gion w a bs < uj < uj c they are meta-stable. For uj > uj c we 
have unstable Q-balls. 

In the limit uj — > w max = m one observes M — > mQ, 
sec Fig. |U This means the unstable Q-balls dissociate 
into Q-clouds, i.e. into a gas of free quanta [28|. 
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FIG. 5: Ground state properties of Q-balls as functions of u. (a) charge Q, Eq. ((6}. (b) mass M, denned before Eq. (|14|l , (c) 
constant di, Eq. (|18|) . (d) "surface tension" 7, Eq. (|34[) . The mean square radii of (e) charge and (f) energy densities, Eq. ((36} , 
and (h) shear forces, Eq. ((34}. (g) position _Ro of the zero of the pressure. The values of densities in the centers of Q-balls for 
(i) charge density, (j) energy density, (k) pressure, Eqs. ((61 1141 [17} . (1) The "surface energy" defined in (|27(l . The special value 
^abs ~ 1.245 (uj c ~ 1.38) is marked by a square (star). The vertical lines indicate the limits oj m i n w 0.447 and oj max « 1.483. 



D. Ground state properties 

In this Section we study "global" Q-ball properties, ap- 
propriate integrals of the "local" densities from Sec. MI Bl 
The numerical results are shown in Fig. [5] which is orga- 
nized as follows. The columns show as functions of u> 
quantities associated with (from left to right) the dis- 
tributions of charge, energy, pressure and shear forces. 
Values of w a bs, w c , w m m/max are indicated in all plots. 

Figs. (S^r-d show Q, M, d\, 7. At ui = uj c charge Q(uj) 
and mass M{uS) exhibit global minima, see ((T2|) . while in 
the vicinity of uj ~ uj c the "surface tension" 7(0;) exhibits 
the largest curvature, and —di(ui) a global minimum. 

Figs. (S^-h show the different length scales of Q-balls: 
square roots of the mean square radii of the charge and 
energy densities and shear forces, Eqs. ((Ml 155} . and Ro 
which is where p(r) changes sign. The behavior is quali- 
tatively similar: all radii increase as uj — > oj m i n / max , and 
have global minima around uj w w a bs (for the minimum 
of (r 2 ) 1 / 2 we cannot exclude that it is, within numerical 
accuracy, exactly at w a b s )- 

Figs. [5J-k shows the charge density, energy density, 
and pressure at the center of the Q-balls as functions of 
uj. These quantities exhibit maxima around u; a bs- 

Since s(0) vanishes, it would make no sense to show 
this quantity in analogy with Figs. 0-k. Instead, in 
Fig. [5j we show the "surface energy," Eq. (f27|) . as func- 
tion of ui. As uj increases from cj m i n to w max the sur- 
face energy decreases monotonically, which is compatible 
with the view that with increasing uj the Q-ball becomes 
a more and more diffuse object, see Sec. MI Bl such that 
the role of a "surface energy" become less and less impor- 



tant. E SUI f(uj) changes the curvature at the point uj = uj c 
within numerical accuracy. 

Some of the quantities vary strongly with w, for in- 
stance di extends over 12 orders of magnitude. This 
is not surprising since we compare Q-balls with dif- 
ferent masses and sizes. For each quantity one could 
find "natural units" in order to make (from this point 
of view) the comparisons quantitatively more meaning- 
ful. The dimensionless constant d\ has the natural units 
(mass x length) 2 . To see this notice that it can be ob- 
tained from di = \M Jd 3 xr 2 s(r), and Jd 3 xs(r) has di- 
mension mass, since s(r) and loo(r) have the same di- 
mensions. Thus, one choice of "natural units" to mea- 
sure di could be M 2 (r 2 ) with i = s,E. In Fig. |6] we 
see that in these units the constant d\ varies much more 
moderately. In fact, for Q-balls of all uj we find that 
< (-di)/(M 2 (r 2 )) < 0.3 holds. 




FIG. 6: (— di) in units of M 2 (rj) (i — s, E) as function of oj 2 . 
These units take into account the true dimensionality of d\. 



8 



IV. THE LIMIT w 



these radii to coincide for lu — ¥ w m ; n , i.e. 



In this section we discuss Q-ball properties in the limit 
lu — > cj m in, the so-called "thin wall" limit. In this limit 
the solutions describe objects of increasing size R with 
uniform charge distribution for r < R, which drops to 
zero over a narrow transition region ("thin wall") [25| . 
In some sense the Q-balls resemble liquid drops. 

In a liquid drop of the size R the pressure distribution 
is p(r) = po &(R — r) — ^p RS(R — r) where po denotes 
the constant pressure inside the drop. This p(r) satisfies 
the stability condition (|22|) . The shear forces are given 
by s(r) = jS(R — r), and the differential equation ([20|) 
leads to the Kelvin relation 7 = ^poR. 

In the following we will "test" the predictions from 
the liquid drop picture using our numerical results, and 
derive analytically relations valid in the limit lu — > Lu m { n . 

In Fig. [2] we have seen that the solutions and densities 
approach the expected liquid drop shapes, see Sec. IIIIBI 
Let us highlight here the shear force distribution. Fig. [7] 
shows s(r)R /^/ as function of r/R in the "edge region" 
for selected lo close to w m i n . The curves are scaled such 
that the areas under the graphs are normalized to unity. 
Clearly, as lu approaches u; m ; n the shear force distribu- 
tions peak more and more strongly in a narrow region 
concentrated around r/Ro sa 1. The "edge region" makes 
makes 5 % and less of the size of the Q-ball, as expected. 

That the Q-ball size diverges as lo — > Lu m - m is apparent 
from Figs. [S^-h. Our first quantitative expectation is 
that both the radius i?o describing the position of the 
zero of the pressure and (r 2 ) 1 ^ 2 characterize equally well 
the position of the "edge" of the Q-ball. Hence we expect 
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FIG. 7: s(r)J?o/7 as function of r/Ro in the "edge region," 
for selected values of lu in the range 0.216 < lu 2 < 0.25. The 
curves are scaled such that the areas under the graphs are nor- 
malized to unity. The figure shows that with lu 2 approaching 
^min = 0-2 the shear forces approach their liquid drop limit 
s(r) — 7 5(r — Ro) where 7 denotes the surface tension, and 
7?o the position at which the pressure vanishes. 



The numerical results in Fig. [5^, support Eq. (|4"2")l . 

In the following we choose (r 2 ) 1 ^ 2 as a reference length 
scale for the size of the Q-ball in the liquid drop limit, 
and define the "surface" and "volume" of a Q-ball as 



A, 



4-7T (r: 



V. 



4-7T 

~3~ 



r 2\3/2 



(43) 



(a) 










R§/< r 2> 






<iQ>/<r^> - 






<r|>/<r^> - 
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(b) 



The charge distribution becomes p c h{f) = po&(R — r) in 
the liquid drop limit, yielding (r 2 h ) ~ \R 2 for the mean 
square charge radius in Eq. (|36j) . The situation is analog 
for (r%), although Too(r) has a i5-function-type "spike" 
at r = R due to surface energy, as can be seen in Fig. [2};. 
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FIG. 8: Q-ball properties as functions of lu 2 in the limit 
lj -> tJmin. (a) the ratios R^/(r 2 a ), {r 2 Q )/(r 2 ), (r|)/(rf). 
(b) M/{T 00 {Q)V S ) and Q/(p ch (0)V s ). (c) 2<y/(p(0)(r 2 ) 1/2 ) and 
{-dx)/{M{M-LuQ){r 2 s )). (d) p ch (0), r oo (0),p(0). The solid, 
dashed, dashed-dotted lines show our numerical results, which 
approach the predicted limits marked by symbols. In Fig. [HJi 
the thin lines are the analytic results derived from Eq. (|50l) . 
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But the surface energy is proportional to R 2 , while the 
contribution of the constant bulk matter density inside 
the drop is proportional to R 3 , so the influence of the 
spike can be neglected for a large drop. Hence, we expect 



lim 



3 

5' 



lim 



(44) 



and the numerical results in Fig. [8^1 support this. For the 
above discussed reasons we can furthermore expect 



lim — 

U-yu^ T 00 (0) Vg 



lim 



Q 



Pch(0 )V a 



= 1, 



= 1 



(45) 



which is also confirmed, see Fig. \E]p. 

Surface and surface tension 7 are abstract notions for 
arbitrary Q-balls which are defined through Eq. (f3~4"|). 
One way to check the usefulness of those definitions pro- 
vides the Kelvin relation, which implies the expectation 



lim 



2- 



MW 2 



= 1 



(46) 



Notice that because of our definitions (|3~4"1 143|) and (pOf , 

we always have the relation 



2-) 



2(M -loQ) 



Po{r 



2U/2 



PO V s 



(47) 



drop 



From (|19|) we obtain for a liquid drop d\ 
Inserting here the expression for 7 from (|47[) yields 



-4fM 7 i? 4 . 



(-l)di 
«-£L M(M-cjQ)(r 2 ) 



(48) 



Fig. confirms both relations (|4"6")l and 

Next let us focus on the center properties of Q-balls. 
For uj — > w m i n the limiting value of the field 0(r) at r = 
assumes the value </> 2 onst = B/ (2C), see (f4*Tj) and App.lA"! 
For our potential and choice of parameters this means 

lim p ch (0) = w min 4> 2 const = V02, 

W-^Wmin 

lim Too(0) = w 2 lin (j) 2 const = 0.2 , 

UJ—^-UJmin 

lim p(0) =0, (49) 



The result for (^const(w) follows from Eq. 
the plus sign, and can be written as 



6 2 

const 



using 



(50) 



Fig. [5U shows that ([50)1 provides excellent approxima- 
tions for the exact w-dependencies of p c h(0), Ibo(0), p(0) 
from up to a; 2 < 1. The reason for that can be seen in 
Fig. [1] In our potential, the "particles" (in the "particle 
motion" picture of Sec. MI Al) are released very close to 
the respective maxima of U c s, and this is what Eq. (|50[) 
actually describes, see App. [X] 

This brings us to another test of the liquid drop pic- 
ture: asw-> w m i n we expect the "edge" of the Q-ball to 
become more and more well-defined. We have seen this 
in Fig. but this observation can be made quantitative 
by defining the thickness of the edge region as 

(Ar 2 ) 2 = (((r 2 - ((r 2 ))) 2 )) = ((r 4 )) - ((r 2 )) 2 > (51) 

where we introduced, c.f. App. [Cl 

f^°dr r n s(r) 



«r n » 



J °°dr s(r) 



(52) 



With this definition we can formulate the expectation 
that in the "thin wall limit" the relative size of the "edge 
region" vanishes 



Ar 2 
hm 7~2T 



0. 



(53) 



This is supported by the numerical results, see Fig. [9] 

Finally we turn our attention to the behavior of inte- 
grated quantities for uj — > w m i n . In this limit Q, M, d\ 
and the radii diverge, as shown in Fig. [5l although cer- 
tain ratios of these quantities remain finite, and follow 
the predictions from the liquid drop picture, see the dis- 
cussion above. The key to understand quantitatively the 
behavior of these quantities in the limit uj — > cj min is the 
surface tension 7, the only "integrated" Q-ball property 
which remains finite in this limit. 



Ar s /<r_> 



which is supported by our results in Fig. [5Ji. The result 
for p(0) is derived alternatively in App. [Bj 

We can go a step further and derive predictions from 
the liquid drop picture for the densities in (|49| also for 
ui 7^ Wrnin. This can be done because the probably most 
important features for uj > w m i n are the finite size of 
the Q-ball, and its diffuse "edge." But these features 
become important "far away" from the Q-ball center. So 
one would expect the liquid drop approach to give useful 
approximations for p c h(0), Tqo(0), and p(0) not only for 
ui = uj min but also for uj in some vicinity of o; min . 
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FIG. 9: The ratio /±r1/{r 2 s ) characterizing the relative size 
of the "wall thickness," as function of uj 2 (solid line). In the 
limit U) — > CJ m in ("thin- wall limit") Ar 2 /{r 2 } -> (marked by 
the symbol). The numerical results support this expectation. 
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FIG. 10: Q-ball properties as functions of u> 2 plotted in the 
form X(lu 2 — uj 2 nin ) N . The respective quantities X, and their 
scaling powers N, written as the pairs (X, N), are as follows: 
(a) ( 7 ,0), (b) «r?},2), (c) (Q,3), (d) (AT, 3), (e) (di, 7). The 
triangles mark the analytical predictions from Eqs. (1571 1591 
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The surface energy i? sur f diverges because it is propor- 
tional to the surface area, which diverges for u) — > oj m 'm- 
Taking out carefully these divergences allows one to de- 
fine the surface tension, cf. Eq. (2.19) in (25[, as 



lim 7 = lim 



2U, 



(54) 



where U = V(<fi) — \ui 2 <\> 2 . Let us define 



With the substitution 



> 0. 



we obtain 



(55) 




For £ m j n 7^ the integrand is complex. Recalling that 

for e min we have </>o(e min ) -> </>const = B /( 2C )> see 
App. fffl we obtain, for our parameters, 



lim 



1 



yc_ (b_ y 



i 



2^2 \2CJ 2^2' 



(57) 



which agrees with the numerical results, see Fig. [TOk . 

Next we want to determine the behavior of the mean 
square radius (r 2 ) in this limit. From Eq. (|50p we obtain 
for p(0), Eq. dTTJ), the behavior 



p(0) 



1 



2 i2 



^(00) 



4C* mm 



+ G(4 lin ), (58) 



which we have seen in Fig. [HJi and derived alternatively 
in App. [Bj Combining this result with (|4l?I |57|) yields for 
our potential 



D2 



(59) 



which is supported by the numerical results in Fig. 110b . 
Analogously we obtain 



7T B 4 
7T 



UJr, 



lim 



r 14 

■■min 



3V2 c 5 / 2 min ' 

144 C 6 Wmin ' 



(60) 
(61) 
(62) 



The numerical results in Figs. [T0b-e fully support these 
conclusions. We see that among the quantities in (f57H6"2"j) 
di has the most rapid rise for u — > w m ; n , which explains 
the observations in Fig. [5l Combining (|59l [BTl [62]) yields 



di 



M 2 (r 2 



1 



4 ui 



(63) 



where the dots indicate higher order terms. Eq. (|63|) 
explains the observation made in Fig. HI namely that d\ 
measured in its "natural units" vanishes in this limit. 

The liquid drop analogy was very successful. One could 
be tempted to drive it further than we did it here, e.g., by 
giving the drop also a uniform charge distribution. The 
resulting repulsive forces ensure stability, and a virial the- 
orem analog to (f29|) can be derived. But the microscopic 
details of the stabilizing dynamics are different from a 
Q-ball, and we will not pursue this analogy further. 



11 



V. THE LIMIT uj ^ max 

In this section we discuss the properties of Q-balls for 
lu — > u> max . For certain potentials one obtains small and 
stable Q-balls in this so-called "thick- wall" limit [30j . 
But in our potential for uj 2 > u> 2 pa 1.9 the Q-balls are 
unstable. For instance, the solution to 2 sa 2.192 with 
Q = 42 and M pa 62.4 can decay into 2 absolutely sta- 
ble Q-balls corresponding to uj 2 pa 1.223 with Q = 21 
and M pa 28.5, or into 3 absolutely stable Q-balls corre- 
sponding to uj 2 pa 1.466 with Q = 14 and M pa 20.4. 1 

Finally, as uj — > w max in our potential, the solutions 
get more and more spread out, and approach from above 
M — > mQ where m = w m ax is the mass of the quanta, see 
Fig. HI This means that the unstable Q-balls dissociate 
into a gas of free quanta, a "Q-cloud" [281 ) . 

The aim of this section is to study analytically how 
Q-ball properties behave for uj — > uj max . The key for 
that is the large-r asymptotics of 4>(r) derived in Eq. (flU)) 
which shows that as long as uj 2 < w max the solutions 4>(r) 
decay at large r fast enough to ensure the convergence 
of the integrals appearing in M, Q, or other properties. 
Of course, the existence condition ([7]) requires uj to be 
always smaller than w m ax- But we may study the scaling 
of Q-ball properties as ui approaches w max from below. 

Let us define 



»(r) 



- uj 2 > 0, 



_ (64) 

which is the leading term in the large-r asymptotics of 
4>(r) in ((TO]) . We consider first the charge Q in Eq. (|6|) 



Q = A-kui J drr 4> (r) 
o 



. + 4ttw / drr 2 (f>l symp {r) 



■ J dx cxp(— 2x) , 



(65) 



where in the second step we split the integral into an in- 
ner (indicated by the three dots) and an outer part. It 
is understood that this decomposition is done at a suffi- 
ciently large radius R such that (j)(r) can be well approxi- 
mated by its asymptotic form for r > R. In the third 
step in (|65[) we made the substitution r — > x = ff max r. 

From Eq. (|65p we see what happens as e max decreases. 
The inner part indicated by the three dots in (|65[) gives a 
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FIG. 11: The Q-ball properties di, M, Q, Ro, n = (r?) 



2x1/2 



(i = E, Q, s), E sul{ , 7 as functions of e max = V^Lx - uj 2 . 
The region in the plot covers the range 1.8 < uj 2 < 2.195. 
The small-e max scaling of di, M, Q, (rg) 1 ^ 2 and (r|} 1/ ' 2 was 
predicted analytically in Eqs. 



finite contribution to Q, but the outer contribution scales 
like l/e m ax- Thus, we expect that with decreasing e max 
the product e m axQ — > const. This method though does 
not allow us to determine the value of the constant. For 
that a more careful analysis is needed, which we will re- 
port elsewhere [46| . But in this way we correctly predict 
that Q oc l/e max at small e max , which is fully supported 
by the numerical results, see Fig. [TT] 

Applying this method to other quantities we obtain 
the results summarized below (all constants are positive 
and different in each case) 



lime max Q = 


const, 


(66) 


lime max il/ = 


const, 


(67) 


!im e max d l = 


—const, 


(68) 


Km £ max ( r k) = 

E— >0 


const, k = Q, E, s, 


(69) 


liln , Smax-^surf = 


const , 


(70) 


lin l£max7 = 


const . 


(71) 



Here we content ourselves to state that the decays are possible 
energetically, but we are not concerned about their dynamics. 
Notice also that in these examples integer charges were chosen. 
But in general the charge Q is not quantized. Also "asymmetric" 
decays into Q-balls of different charges are possible, but then less 
energy is released. 



The predictions (|66H71|) are fully supported by the nu- 
merical results as shown in Fig. [TTJ 

Notice that the results (foWTTj) for E smil 7, (r 2 ) are 
numerical observations, because our method cannot be 
applied to quantities vanishing with e max — > 0. In fact, 
for instance the scaling of the outer contribution in the 
integral (|2"7|) defining E SUT { does imply Eq. (|7D|) . But our 
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rough method would gcncrically suggest that the contri- 
bution of the inner region scales like e^ ax and dominates. 
A more careful analysis is needed to prove the prediction 
(f70|) , see |46| . The same reservations apply to the scaling 
behavior of 7 and (r^), which are both connected to -Esurf 
via Eq. (|53|>. 

We observe numerically that the position R where 
p(r) changes sign scales in the same way as the square 
roots of the mean square radii in (|69p . Thus, indepen- 
dently of whether we measure it in terms of Rq or the 
(r 2 ) 1 / 2 , the size of the solutions grows with e max — > 0. 

The constant d\ diverges as l/e^ax with decreasing 
£max, sec Fig. 1111 However, when measured in its natural 
units it actually goes to zero as di/(M 2 (r 2 )) oc e^ax as 
we have seen previously in Fig. [5] 

To summarize, as e max —¥ mass, charge and size of 
the Q-balls diverge as l/e max , see ([551 1571 . Thus, the 
mean charge and energy densities, which are proportional 
to Q/(size) 3 and M/(size) 3 , vanish like e^ax- Fig.g]has 
shown that the Q-balls are unstable, and their (positive) 
binding energy M — mQ approaches zero (from above) 
as £ max — > 0. Hence, in this limit we obtain a dilute gas 
of free Q-quanta as discussed in i28j . 



of the energy functional. But the other way round, the 
pressure is ultimately related to the issue of stability by 
Eq. (|22|) . which we have shown to be equivalent to the 
virial theorem (|29|) . A proof that d\ < on the basis of 
p(r) is therefore in general on a much more solid ground. 



A. Arguments based on s(r) and inequalities 

In this section we will show that d\ < using argu- 
ments based on the shear force distribution. The argu- 
ment is trivial and makes use of the observation that 
manifestly s(r) = 4>'{r) 2 > Vr. 

In Eq. (|18|) we have seen that d\ is given by (— ^irM) 
times the integral over r 4 s(r) over r from zero to infinity. 
Since s(r) > this immediately implies di < 0. This 
inequality can be improved by recalling that <j>'{r) < 
for < r < 00, see Sec. MI Al Therefore d\ < which 
completes the proof. 

The fact that s(r) > can be further explored to derive 
an inequality showing that d\ must be negative. Using 
flSTJIHl) we have di = -^Af 7 ((r 4 )) and (r 2 ) = ((r 2 )) 
and can rewrite the constant d\ as 



VI. THE SIGN OF di 

In this section we will show in several independent ways 
that d\ is negative. In Sec. IVI Al we will use for that the 
observation that for Q-balls s(r) happens to be positive 
for < r < 00, and in Sec. IVI 51 we will explain why s(r) 
must be positive. In Sec. IVI CI we will prove that d\ < 
using arguments based on p(r) and stability. 

One may wonder why several proofs are needed. In- 
deed, the EMT conservation dictates that s(r) and p(r) 
are connected by the differential equation (|20p . which is 
the origin of the equivalent presentations (|T5ini)|) for d\ in 
terms of s(r) and p(r) [l6l |. and we have explicitly proven 
that our expressions for s(r) and p(r) satisfy (p2"D|) . So, if 
one is able to conclude from s(r) the sign of di, then it 
must be possible to draw the same conclusion also from 
p(r). Therefore, at first glance it may seem sufficient to 
conclude the sign of d± in one way, and below we will sec 
that for Q-balls it is much easier to use s(r) for that. 

However, concluding the sign of d\ from s(r) alone 
bears some danger, because from any "input" s(r) one 
obtains via ([2"0|) a pressure p(r) which automatically 2 
satisfies the stability condition (|22|) [l6| . So one may 
well encounter an approach with s(r) > and conclude 
d\ < without being sure one really deals with a correct 
solution of the equations of motion and a true minimum 



di 



2 The differential equation II2UH allows one to determine p(r) from 
a given input function s(r) only up to an integration constant. 
But the latter can be fixed by demanding that for a well-localized 
finite-energy object p(r) — > as r — > 00. In a similar way one 
can determine s(r) from a given input function p(r). 



M 2 (r 2 




73T ■ (72) 



Notice we implicitly benefited from the fact that s(r) > 0, 
when introducing the averages ({r n }) in ([52]) . Next we 
explore that E aul f = Jd 3 r s(r) > 0, and with w Q > we 
conclude from (J3UJ) that < M - luQ < M. Using the 
latter inequality in (|72[) finally implies that 



< 



di 



< 1 



Arj 

/r-2\ 



(73) 



This proves that (— di) > 0. As a byproduct Eq. ([73^) 
provides also an upper bound on (— di) but in terms of 
Ar 2 /(r 2 ). At this point it is not obvious whether this 
quantity is bound from above, though numerically we 
observe this to be the case in Fig [5] That Ar 2 /(r 2 ) is 
indeed bound from above will be shown in (46j . 

For completeness let us mention the following more 
useful upper bound on the magnitude of (— di). The 
staring point is Eq. (|35|) where we neglect the positive 
quantity uQM (r|) = M d 2 Jd 3 r r 2 cj)(r) 2 . This yields 
the bound 



di 



M 2 (r 2 E ) 



< 



(74) 



which is satisfied by the numerical results, see Fig. [6] 
We checked that this is the strongest inequality one can 
derive involving (r 2 E ) as length scale. The inequality ([74]) 
is atrractive because it provides an upper bound on (— d{) 
in "its natural units" solely in terms of quantities related 
to the energy density T 00 (r). 
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B. Arguments based on the particle interpretation 

In the previous section we explored the observation 
that s(r) happens to be positive for < r < oo for Q- 
balls. Here we will show this must be the case. 

For that we use the particle interpretation picture pBT ] 
discussed in Sec. IIII Al The Newtonian equation (|39|) de- 
scribing the motion of a unit mass particle moving in the 
effective potential U c g = \u) 2 x 2 — V(x) under the fric- 



tion Ff r j 



x(t) follows from the Langrange-function 



L(x,x) and Rayleigh's dissipation function T{x) 1 

L(x, x) = ^ ±{tf - U cS (x), F{£) = j ±{tf , (75) 



according to 



d_fdL\ _9L_ _dF 
dt\dx J dx dx 



(76) 



The physical meaning of Rayleigh's dissipation function 
F{x) is that it describes the rate at which the system 
dissipates its energy E due to the frictional force, namely 



dE d /_ . dL T 
—r- = ~r\ 9x — — L 
dt dt V dx 



= -2T< Vi, 



(77) 



which must be negative because the system dissipates 
energy This means that J-(x) > Vi. 

If we recall that x(t) and t in the particle interpretation 
picture correspond to 4>{r) and r, we instantly see that 
^(x) corresponds to ^s(r). Since -F(x) > this proves 
that the distribution of shear forces s(r) > 0. 



C. Arguments based on the pressure p(r) 

In this section we will prove that d\ is negative, basing 
our arguments on the pressure distribution. 

Let us first demonstrate that for oj satisfying the ex- 
istence condition (JT]) the pressure is positive at small r, 
and negative at large r. In Sec. IIII Al we have proven 
p(0) > 0, and for reasons of continuity p(r) > also in 
some vicinity of the origin. At large-r we derive from 
([TO]) the following asymptotics for the pressure 



p(r) = - (u4 ax -w 2 ) 



2 c 2 
3r 2 



exp 



(78) 

where the dots indicate subleading terms. Clearly, for all 
uj satisfying the existence condition (J7J), the pressure is 
negative at large r. To summarize, we have 



p(r) > for small r, 
p(r) < for large r. 



(79) 



This implies that p{r) must change the sign an odd num- 
ber of times. Of course, p(r) must change sign at least to 
comply with the stability condition (j2"2"j) . From physical 



point of view, we expect p(r) to be positive in the center 
(which implies repulsive forces directed towards outside) 
and negative outside (attractive forces towards inside), 
as we derived in ([79"]) . A stable solution arises when the 
repulsive and attractive forces exactly balance each other 
according to (|22[) . This physically intuitive pattern was 
observed also in soliton models of the nucleon [l6|-[l9[ . 

For a ground state one may expect the pressure dis- 
tribution to change sign only once, see Fig. If we 
assume p(r) to change sign one and only one time, this 
immediately implies that d\ is negative. Fig. [12] illus- 
trates the argument. The left panel of Fig. [T2l visualizes 
the stability condition ([2"2"]l : the shaded areas above and 
below the x-axis are equal and exactly compensate each 
other. Thus, due to the stability condition (|22[) we have 

Ro 



J dr R 2 r 2 p{r) = - J dr Rlr 2 p{r) 




(80) 



which means J^dr r 4 p(r) < 0, and d\ must be negative, 
as can be seen in the right panel of Fig. rT2] 

If we knew p{r) has one zero only, the proof that d\ < 
would be complete here. It is intuitive to assume that the 
pressure distribution of a ground state changes sign only 
once according to (|79[) . However, here we will provide 
a general argument based on stability which is valid not 
only for ground states. 

For that we need the following lemma. For any solution 
of the Q-ball equations of motion we have 



/ 



dr r 2 p(r) > for < R < oo. 



(81) 



To prove ([81]) we make use of the result ([25]) derived in 
Sec. Ill CI and explore the particle interpretation picture 
of the Q-ball equations of motion. For that we notice 
that p(r) + | s(r), the right-hand-side of (1251) . is positive 
V r < oo because 



p( r )+ o s ( r ) 



0'(r) 2 + -w 2 <Kr) 2 - V{4>) . (82) 



-Ekin 



In other words, p(r) + | s(r) corresponds to the total, 
kinetic plus potential, energy of the particle at a given 
time t (with t -f> r). The total energy of the particle must 
be larger than zero V t < oo, because at any finite time 
t the particle still needs to do work against the friction 
in order to arrive, after infinite time, at the origin x = 
(with x <f>). This means that the integral in ([23]) is 
positive and proves Eq. ([8T]) . 

To understand this result intuitively, we remark that 
if the integral ([5T]) could become zero at some finite R, 
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r p( r) /a N 



(a) r p(r) /a N 



(A) 







w 2 =0. 


8 



-A 




(O 2 =0.8 





0.5 




FIG. 12: r 2 p(r) and r 4 p(r) as functions of r for selected U). For 
better comparison, r is given in units of the radius Rq where 
p(r) changes sign, and the normalization factors ajv are such 
that the curves reach unity at their global maxima. The left 
(right) panel shows the integrand of the stability condition 
(the integrand of di, Eq. (|19[) , up the prefactor 57rAf). The 
figure illustrates why di is negative. Integrating the curves 
in the left panel yields zero due to the stability condition. 
Weighting the curves by an additional factor of r 2 and inte- 
grating then yields a negative result for di, see right panel. 



then the fields would stabilize themselves in a subinterval 
r £ [0,R]. Then, setting the fields outside that interval 
to zero, would yield a stable solution with lower mass 
(because Too(r) > also in the omitted region), in con- 
tradiction to the expectation that a given set of initial 
value data leads to a unique minimum of the action. 

We have now all ingredients for the proof that d\ < 
based on pressure and stability arguments, namely ([75)1 
and (|8ip . The proof is as follows. 

Eq. (|79[) means p(r) must change sign an odd number 
TV of times. Let Ri with i = 1, . . . , N denote the radii 
where this happens with 0<i?i<i?2<-..<-Rw<oo. 
Notice that we do not include points where p(r) could 
have zeros without changing sign. 

The stability condition (|22p can then be written as 

Ri R.2 oo 

far r 2 p(r) + Jdr r 2 p(r) + . . . + Jdr r 2 p(r) = . (83) 

Ri R N 

By construction p(r) > in the first, third, . . . integrals, 



and p(r) < in the second, fourth, . . . last integrals. We 
will now replace each of the terms in ([55)1 by a smaller 
term, and show in this way that f °°dr r 4 p(r) < 0. 

Step 1. We consider the first 2 terms in ([53")) . In the 
first (second) term p(r) is positive (negative). Therefore 

Ri Ri 

J dr r 2 p(r) > J dr r A p(r) , 



R 2 



dr r p(r) > 



1 

R 2 



o 

R 2 



dr r p{r) . 



Ri Ri 

Adding up the 2 inequalities in (|84| we obtain 

R.2 R.2 



dr r 2 p(r) > 



1 

Rj 



dr r A p(r) 



(84) 



(85) 



If there is only one change of sign, then we take the limit 
i?2 —¥ oo and recover the situation of Eq. ([50)1 . and our 
proof is complete here. If p(r) changes sign more than 
once, i.e. N > 3, then we have to continue our proof and 
include further contributions in step 2. 

Step 2. Notice that J^ 2 dr r 4 p(r) > for R 2 < oo due 
to (JHTJ). Moreover Rx < R 3 < oo. Therefore 

R.2 R.2 R.2 



1 



1 



dr r 2 p(r) > — ^ / dr r A p{r) > — j / dr r A p(r) . (86) 
R 1 J R 3 



In the next two intervals p(r) > for r € [i?2,-R3], and 
p(r) < for r G [i?3,i?4]. Therefore, in analogy to 



R 3 Ra 
far r 2 p(r) > J dr r 4 p(r) , 
R.2 3 R2 

R.4 i?4 

'drr 2 p(r) > — j dr r 4 p(r) . (87) 

^3 J 

R 3 R 3 
Combining (|86| and the results in (|87|) we obtain 

Ri Ri 



dr r 2 p(r) > 



dr r 4 p(r) 



1 



If p(r) changes sign exactly N = 3 times, then we take in 
(1581 the limit R4 — > 00 and our proof is completed here. 
If p(r) changes sign more often, then we repeat step 2. 

Last step. If p(r) changes the sign N — 2k + 1 times 
(recall that N is odd, and Rn < 00), then we repeat 
successively the 2 nd step fc-times, until we arrive at 



00 00 

/dr r 2 p(r) > — 5- / dr r A p 
Rn J 



(r) 



(89) 



Now, the first integral vanishes due to the stability condi- 
tion (|22|) , and using the definition (fT9|) giving di in terms 
of the pressure, we conclude from ([89]) the desired result 
that d\ must be negative. 
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D. Stability and d\ 



is strictly negative for all Q-ball solutions, 



In the above section we have proven that the stabil- 
ity condition (|22|) and the properties of the pressure at 
small and asymptotically large r unambiguously yield to 
a negative sign of the constant d\. Therefore, clearly if 
we have a stable object then d\ < 0. 

The inverse, however, is not true in general. A negative 
d\ does not imply the object is stable. Fig. IT2"1 illustrates 
this point. The solution in Fig. [12a for ui 2 = 0.8 < w 2 bs 
is absolutely stable. The Q-ball solution in Fig. [T2"b for 
w abs < w 2 = 1.7 < ui 2 is metastable. The solution in 
Fig. IT2"b for ui 2 < uj 2 = 2.0 is unstable. But in all 3 cases 
d\ is negative, as the Figs. [T2H~e illustrate. 

The point is that all solutions for w m in < w < w max 
correspond to minima of the energy functional. There- 
fore all solutions satisfy the stability condition (|2"2"j) . or 
cquivalently the virial theorem (f2"9")l , as well as all condi- 
tions required to prove that they have a negative d\ . But 
only for u) m - m < co < w a bs we deal with global minima of 
the action and absolutely stable Q-balls. 



VII. CONCLUSIONS 

We have presented a study of Q-balls in a scalar field 
theory with U(l) symmetry, and investigated the prop- 
erties of Q-balls as functions of the angular velocity uj. 
While Q-ball stability was studied in literature before 
j37T - l4ll ]. to best of our knowledge this is the first study 
in which this issue is addressed from the point of view of 
the EMT, and the constant d\. All solutions presented 
in this work were exact solutions of the equations of mo- 
tion. Particular focus was put on the behavior of Q-ball 
properties for ui approaching the boundaries w m i n ,max of 
the region in which solutions exist. 

For ui — > u) m in the Q-balls occupy increasingly large 
volumes filled with Q-ball matter of nearly constant den- 
sity [III. We have shown that in this limit the Q-ball 
properties follow the predictions of the liquid drop pic- 
ture. Certain Q-ball properties such as charge Q, mass 
M, mean square radii, and d\ diverge as ui — > Lu m - m . We 
derived analytically the limits for these and other proper- 
ties, which are fully supported by our numerical results. 

In the opposite limit uj — > u> max the solutions become 
unstable and approach the "Q-cloud limit" [28| . Also in 
this limit some properties diverge. We derived analyt- 
ically the scaling behavior of these quantities as uj ap- 
proaches w max . Further results will be reported in [4q [. 
It is remarkable that, among all properties we studied, 
d\ diverges most strongly as ui — > w m m. max . 

The conservation of the EMT implies among others 
the stability (von-Laue-) condition [j, |48| which states 
that the pressure must satisfy / °°dr r 2 p(r) = 0, which 
we have proven explicitly in two independent ways. One 
of the proofs is equivalent to the virial theorem. 

The central result of this work is that the constant d\ 



dx < , 



(90) 



for which we have provided 2 explicit analytical proofs. 
One proof involved the relation of d\ to s(r), and made 
use of the Newtonian particle interpretation of the Q- 
ball equations of motion (25| in which the shear force 
distribution s(r) is related to the Rayleigh dissipation 
function describing the frictional forces. Since the New- 
tonian system dissipates energy due to friction, s(r) must 
be positive for < r < oo. This implies d\ < 0. 

The other proof explored the relation of d\ to p{r), 
and made use of stability arguments. We have shown, 
using the equations of motion, that p(r) is positive in 
the center of the Q-ball (which corresponds to repulsion) 
and negative at large r (which corresponds to attraction). 
This means p(r) must change the sign an odd number of 
times, and we have formulated a general proof valid for 
any Q-ball solution with a pressure with an arbitrary 
number of zeros. We observed that for ground states 
p(r) changes sign only once, but for radial excitations of 
Q-balls one encounters more complex structures (46j . 

The proof of based on p(r) elucidates that for Q- 
balls the negative sign of d\ is a consequence of stability. 
This was conjectured, but could not be p roven rigorously, 
also in other soliton approaches [lg-[l8[ ■ 

The last important insight of our study is that stability 
implies d\ < 0, but the opposite is not true. A negative 
di docs not necessarily mean the object is stable. In 
fact, our proofs hold equally for stable, metastable and 
unstable Q-balls. To guarantee a negative d\ both proofs 
require that we deal with the exact solution of the equa- 
tions of motion, i.e. with a minimum of the action, be it 
global or local. But absolute stability requires that the 
minimum of the action is global. Thus, stability is a suf- 
ficient but not necessary criterion for d\ to be negative. 

An interesting open question is how quantum fluctua- 
tions HU modify the picture of d\ and the proof d\ < 0. 
It would be also interesting to see how d\ is altered for 
Q-balls cou pled to fermionic fields, which allows them to 
'evaporate' [231 ■ The ultimate goal would be to generalize 
the proofs given in this work to quantum field theories, 
and to apply them to the description of hadrons. So far, 
in all theoretical studies diwas always found ne gat ive, for 
pions [3, nucleons ^U^- and nuclei 0, [H H2 . Also 
lattice QCD calculations yield a negative quark contribu- 
tion to the d\ of nucleon, though the gluon contribution 
and hence the total d\ are not yet known [l4[. First ex- 
perimental results are compatible with d\ being neg ative 
[50j but this observation is not yet conclusive [5l|, and 
future data will provide further insights (52j . 
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Appendix A: Stationary solutions 



Appendix B: Bounds on the pressure 



In this Appendix we discuss stationary solutions of the 
equations of motion (jU [5]) of the type <f>(r) = CO nst Vr. 
Though they do not obey the boundary condition ([5]) at 
infinity, these solutions are nevertheless of interest. 

If <fi(r) = C onst the boundary conditions at r = in 
([5]) hold, and the equation of motion (0} becomes 



The pressure at the origin is just p(Q) — U g{4>o) an d 
this is positive, Eq. (|40|. because for any regular solution 
the effective potential of the particle at the starting point 
must be positive. In this Appendix we will show that 
the pressure is also bound from above. From ([5j I17[) we 
obtain 



W 2 0const - V'((f> CO n 

uj 2 -2A + AB<t>, 



(Al) 



const 



6C# 



const 



The trivial solution </> CO nst = 
Further (for uj m i n < uj < w max 



describes the vacuum, 
always real) solutions are 



b 2 

const 



B 
3C 



± 



B 2 
9C 2 



'■-2A 
6C 



1/2 



(A2) 



P(0) = j w >o 



1 2j?fi 2 ^0) 
2^2 



W 2 </> 2 



(Bl) 



Using © we notice that V <fi(r) (including <j>(r) at r = 0) 



> min 



2U(0) 



(B2) 



Inserting (|B2[) in (|B1[) and including also the lower bound 
from Eq. P0|) we obtain 



Because £ in ([T]) is symmetric under </>—»—(/> it is suffi- 
cient to focus on the non-negative solutions. 

In the particle interpretation picture the two positive 
stationary solutions in (|A2[) have the following meaning. 
The solution with the minus-sign in (| A2|) corresponds to 
the situation that the particle is at t = precisely in 
the local minimum of the effective potential in Fig. [T] and 
will stay there forever. The solution with the plus-sign 
in (|A2[) corresponds to the (not stable) situation that the 
particle is at t = precisely in the global maximum of the 
effective potential in Fig. [1] and will stay there forever. 

Asw-> uj min we obtain </> 2 onst = B/(6C) from (|A2D . 
which corresponds to the (not interesting for us) situation 
where the particle stays forever in the minimum of U e g, 
and the (much more interesting) solution 



b 2 

const 



B 
2C 



= 1 as uj — > w n 



(A3) 



(B3) 



0<p(0)<i(w 2 -< in )^ 



An important application of the upper bound in (|B3|) is 
that it allows to verify independently the result for p(0) 
in the limit uj — > w m i n quoted in (|49p . We remark that 
the upper bound in (|B3[) is for w m i n < uj < oj max always a 
real inequality and saturated only in the limit uj — > oj m i n . 



Appendix C: Generating functional for ((r™)) 

In Scc. lIVI wc defined the ((?'")) in ([52^1 . which allowed 
us to express compactly (r 2 ) = ((r 2 )) and the measure 
for the wall width (Ar 2 ) 2 = ((r 4 )) - ((r 2 )) 2 . Another 
interesting application, if we continue to negative n, is 
p(0) = 27 ((r -1 )). This allows us to express the result 
obtained in (l46l) as 



which corresponds to the situation with the particle 
placed at the maximum of £7 e fr(<^) which for uj approach- 
ing w m i n from above is just above zero. This situation is 
of interest, because as uj — > uj min the particle has to be 
placed very close to this maximum of U e s(<j)), and "wait" 
there long enough such that its small initial potential en- 
ergy U c ff(4>) is sufficient to overcome the time-dependent 
friction which decreases with time 25 1. 

As uj — > cj max we obtain from (|A2[) the solutions 
2 onst = and 2 onst = (2S)/(3C). When uj = uj max 
the effective potential does not dip below zero at all, i.e. 
it is not possible to release the particle from any <fio > 
so it would stop in the origin [25j. The only solution is 
'/'const = 0- However, from (|A2[) that solution develops in 
the limit uj — > w max from the minimum of U e fi which is 
below zero. But for any regular solution with uj < cj max 
the potential at the starting point U c s((/)q) > 0. There- 
fore, the stationary solution </> 2 onst = obtained here is 
not continuously connected to the limiting value for (j)Q 
stated in Eq. ([Hjl. 



lim ((r- 1 )) ((r 2 )) 1/2 = 1. 



(CI) 



The result in Eq. (|C1|) can be understood and interpreted 
by recalling that s(r) — > 7 6(r—R) in the liquid drop limit 
which is equivalent to uj — > w m i n , see Sec. IIVI 

The positivity of the shear forces, which was crucial in 
Sees . I VI Al and I VI Bl allows us to introduce the functional 



(C2) 



F(X) = J dr s(r) exp(-Ar 2 ) 


which is a generating functional for ((r n )} for even 



00 (_-\\n 

F(A) = F(0)^^((r 2 "))A" (C3) 



The surface tension is just 7 = F (0), the surface energy 

3 



is E suxi = -4ttF'(0) and d x = -4f M F"(0). The mean 
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square radius of the shear forces, and the measure of the 
wall width Ar 2 can be expressed as 



(r 2 s ) 
(Ar 2 ) 2 



d_ 

dX 

dX 2 



log F(X) 
log F(X) 



A=0 



(C4) 
(C5) 



We remark that each integral over s(r) can be traded 



for an integral over p(r) by exploring the differential 
equation ((20)) . For instance, 7 = f / °°dr p(r) and 



(r 2 s ) 



J^dr r 2 p{r) log r 
J °°dr p{r) 



(C6) 



These and further relations were derived in [16j and can 
be used as cross checks for the numerics. 
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